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We present the general form of potentials with two given energy levels 
Ei, E2 and find corresponding wave functions. These entities are expressed in 
terms of one function £(x) and one parameter AE = E%-E\. We show how the 
quantum numbers of both levels depend on properties of the function 
Our approach does not need resorting to the technique of supersymmetric 
(SUSY) quantum mechanics but generates the expression for the superpoten- 



The potentials whose spectrum can be found exactly are very rare in quantum mechanics. 
Meanwhile, the condition of exact solvability can be weakened: one may demand that only 
for a finite part of the spectrum eigenstates and eigenvalues be found explicitly or from a 
finite algebraic equation. This opens two different possibilities. First, there exist so-called 
quasi-exactly solvable (QES) systems, whose Hamiltonian can be expressed in terms of 
the generators of the algebra having finite-dimensional representation (for one-dimensional 
potentials the relevant algebra is SI2, the corresponding generators having the meaning of 
the effective spin operators) [0-0]. In so doing, the dimension of the finite subspace of 




tial automatically. 



PACS numbers: 03.65.-w, 03.65.Ge 



I. INTRODUCTION 




the whole Hilbert space is determined by the value of the effective spin that usually enters 
the QES potential as a parameter. Second, instead of relating the dimension of the finite 
subspace to an underlying structure of a Lie algebra representation, one may fix the number 
of known levels "by hands". In the simplest case this number is equal to two, so we deal 
with two-dimensional subspace. Although such a procedure makes the underlying algebraic 
structure more poor, it extends considerably the set of potentials with the known part of 
the spectrum. 

Physical motivation for interest in potentials with two known energy levels stems from 
the fact that a two-level system represents a very wide class of models often used in solid 
state and nuclear physics and quantum optics. Let us mention here only few examples: the 
Dicke model of interaction between atoms and radiation ||, Lipkin-Meshkov-Glick model of 
interacting nucleons || , the phenomenon of macroscopic quantum tunnelling ]7) . We would 
like to stress that it is just the potential description of systems with a finite number of energy 
levels that enabled one to give clear and simple explanation of the phenomenon of spin tun- 
nelling ||. Therefore, finding potentials, that correspond to a fixed numbers of eigenstates, 
was an important step in calculation of tunneling rates in ferro- and superparamagnets. 

Meanwhile, there is also the inner motivation that stems from quantum mechanics as 
such. From general viewpoint, recovering potentials from a known set of eigenvalues is noth- 
ing else than the reduced variant of the inverse scattering problem. As is well known, using 
Darboux transformation, one can get many-soliton solutions of the Schrodinger equation 
with N energy levels, fixed in advance. Understanding, how the truncation of the scattering 
data modifies the structure of the theory, could gain further insight into the inverse scat- 
tering approach. The first necessary step here is to find the full solution of the problem for 
N=2. 

If N=l (only one level is fixed), it follows from the Schrodinger equation that the potential 
is U — E + ip" /if), where E is the value of energy, if) is a wave function. Choosing any ip(x) 
having no zeros at the real axis, we obtain immediately the corresponding potential U(x), 
regular on the real axis. We would like to stress, however, that, in contrast to the N=l 
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case, when the solution of the problem is straightforward, already for N=2 resolving this 
problem needed the elaboration of different approaches discussed in literature. The existence 
of exact solutions with two levels for power-like potentials was indicated in ||, [Kj. The 
rather powerful technique based on supersymmetric (SUSY) quantum mechanics (see the 
review |nj) was suggested in jnfl , [Tj|. It enables one to generate the potentials with known 



ground and first excited states. The aim of the present paper is to suggest a general approach 
to the potentials with two known levels valid for any n-th excited states. The corresponding 
method and results turn out to be surprisingly simple and do not require sophisticated 
technique (for instance, such as SUSY quantum mechanics). 

II. BASIC EQUATIONS 

Consider the Schrodinger equation with the Hamiltonian H = — + U(x). Let ib\ and 
ib 2 be wave functions obeying the Schrodinger equation: 

= E^x, (1) 
Hib 2 = E 2 i\) 2 . (2) 



Then it follows from ([!]), @ that 



ib" 

U = E 1 + ^, (3) 



E 2 

1p2 ^1 



(4) 



Let, by definition, 



ib 2 = tib x . (5) 

Then we have for ib 2 from @: 

where AE = E 2 — E%. By substitution of (^) and (|5|) to ([|), we obtain three equivalent 
forms for the potential: 



U = E 1 - ^ + - + + i(AE) 2 (i) 2 , (7) 

^ = ^ " ^(1 " 2^) + \(&E) 2 (p 2 - \[t]*, (8) 
U = E x + x' 2 - X", (9) 



where [£\ x = ^ — f (§^) 2 * s Schwarzian derivative (see, e.g., Ch. 2.7 of Ref. [fL4"| ). The wave 
functions of the states under discussions are 

V>i = e" x , ^2 = e"*£. (10) 

Eq. (0) gives us the general formula for the potential with two given energy levels. It 
is expressed directly in terms of their values Ei, E 2 as parameters and one function £(x), 
corresponding wave functions are given by (||), (pi]) and expressed in terms of the same 
quantities. It is worth noting that the function £ (x) does not enter the set of known data - 
rather, the freedom in its choice reflects the fact that for two given eigenvalues there exists an 
infinite number of potentials having two fixed eigenvalues. The structure of these potentials 
is not arbitrary but is governed by the form of £(x) according to (7). 

Eqs. (0), (|5|) and fliTf ) constitute the main result of this paper. It is worth stressing that 



the derivation of eq. (0) is very simple, direct and does not need sophisticated technique, 
such as SUSY machinery. On the other hand, the potential in terms of the function x' has the 
form typical for SUSY quantum mechanics, automatically. In so doing, x' plays the role 
of a superpotential. As is well known (see, e.g., |llj), one- dimensional quantum mechanics 



can always be formulated in a SUSY way. However, given a potential, the superpotential 
cannot be found explicitly for a generic model. Meanwhile, in our case we found not only 
the potential but the explicit expression for the superpotential (|j) as well. 

It is worth stressing that the derivation of (0) - ([|) relies strongly on the successful choice 
of the function ^(x) that parametrizes the family of solutions. The fact that, for given 
Ei, E%, the ratio of two eigenfunctions determines the potential completely generalizes the 
observation made in Refs. ||12||, JTB|| for the particular case when the eigenfunctions under 



consideration refer to the ground and first excited states. 



The formalism elaborated above for the one-dimensional Schrodinger equation can be 
also applied to the three-dimensional one for a particle moving in a spherically-symmetrical 
potential U (r) . After the separation of variables, the effective potential entering the radial 
part of the Schrodinger equation, is equal to U e f — U+ . Then, repeating calculations 
step by step, we obtain 

u-u^ + x* e A V + e "i \ AEe I A ~ 2 ^i + 1) (u) 
u ~ u + A 4^-^|7(- + y)- A ^^- + w > C 11 ) 

\=(l 2 -h)(l + h + l 2 ) (12) 

and is expressed in terms of E\, E 2 and £ by the same formulas (0) - @ as in the 
one-dimensional case. 

It is worth noting that now a new interesting possibility can arise that is absent in the 
one-dimensional case: AE = 0. It becomes possible due to the fact that two quantum states 
can refer to different effective potentials (Zi ^ l 2 ): we are faced with degeneracy with respect 
to the angular momentum. Then the potential acquires the form 

U - El ~ 2* + A 4^ ~ t^' { ~t + J 1 + 2? ' (13) 

We will not discuss the three-dimensional case further and will concentrate on the one- 
dimensional one. 



III. GENERAL PROPERTIES AND CLASSIFICATION OF STATES 

The potential U = U(E 1 , E 2 ,C,) possesses the symmetries that follow directly from (0): 

U(E 1 ,E 2 ,aO = U(E 1 ,E 2 ,0, (14) 
U(E 1 ,E 2 ,£) = U(E 2 ,E 1 ,r 1 )- 

Throughout the paper we assume that the potential U(x) is regular everywhere, except, 
perhaps, infinity. Then all zeros and poles of the function £(x) are simple - otherwise the 
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potential U would become singular and the wave function %p2 would cease to be normalizable. 
If the function £ has a pole at x = xi, £ ~ A(x — xi)^ 1 , one gets from (j^) that x' ~ 
—(x — xi)^ 1 , so ipi(xi) = 0, ip2(xi) = const ^ 0. Therefore, every zero of £ generates a node 
of the wave function ip2 and every pole of £ generates a node of ipi. 

The set of possible nodes of wave functions depends also on the behavior of the function 
x(k) m the vicinity of zeros of the function £'(x) due to possible zeros of the factor exp(— x) 
in (|10|). Let £'(#o) = 0. Then, according to if x — > xo, x' ~ -B/(x — xo), where £> = 
^ "^g^fo)^ ' anc ^ ^ ne P°t en ti a l contains the term that behaves like B(B + l)(x — xo)~ 2 . The 
regularity of the potential entails B — or B — — 1. Consider these two cases separately. 

Let, first, 5 = 0. Now the condition 

(e" + A£?O|x=x O = (15) 

must hold. In so doing, the function is regular in the vicinity of xo due to the 
condition (|I5D and the factor exp(— x) cannot vanish. 
Consider the case B — — 1. Now we have 

(3£" + A££) 1^ = (16) 

Then » — (x — xo) -1 and the functions and ^2 share the common node at x = xq 
due to the factor exp(— x), as it follows from fllPl) . (For example, if the potential is even, 
U(—x) = U(x), all wave functions of odd states vanish at x — 0.) 

As a result, we arrive at the conclusion that, if (i) £(x) has n\ poles and n,2 zeros, (ii) £'(x) 
has m(°) zeros such that (|15"D is satisfied (B = 0) and m* - ) zeros such that ([16j) is satisfied 
(B = —1), the function ipi(x) describes the state with the number of nodes N\ = n\ -\-mS~\ 
while ?p2(x) corresponds to the state with the number of nodes N 2 = n 2 + rm~\ Therefore, 
the quantum number that label states is equal to N\ for ipi and N 2 for ip 2 (^1,2 = 0, 1, 2...). 

It follows directly from the definition (^]): if has simple zeros at X{ and ip2 has simple 
zeros at Xk with x, 7^ Xk, the function £ has poles at x = Xj and zeros at x = x^. However, if 
some Xj = Xk, corresponding zeros of both functions compensate each other and this results 



in the fact that, if some coefficients B = —1, the state labels are not determined completely 
by the numbers ni, n 2 . 

Let us have two fixed energy levels E\, E 2 (E 2 > E\) and the function £(x) such that 
N 2 > Ni. Then, it is the potential U(Ei, E 2 ,£) for which the level E\ belongs to the iVi-th 
state and E 2 corresponds to the A^ 2 -th state. If Ni > N 2 , the relevant potential is U(Ei, E 2 , 
= U(E 2 , Ei t £), the level E\ belongs to the A^-th state, while E 2 corresponds to the 
iV"i-th state. If N\ = N 2 , the function is not suitable for constructing U(x) with two 
different fixed levels. Indeed, in this case we would have had two different wave functions 
with the same number of zeros corresponding to two different levels, in disagreement with 
the oscillation theorem. 



IV. COMPARISON WITH SUSY APPROACH AND TKACHUK'S RESULTS 

Let us introduce the function W + according to 

W + = |, (17) 

where 5 = AE > 0. 

Then, with (^) taken into account, we obtain 

x - = \ (W+ - s W s E±^L, U = E 1 + W>- W, (18) 

fa = e -JdxW^ ^ = £ e -JdxW = W+exp [__ J dx(W + + W-)], 

where, by definition, 

Since fa must be normalizable, sign(iy + (±oo)) = ±1. Let W + have only one zero at x — Xq. 
If we want W_ to be regular at x = x , W' + (x ) = S. 

The formulas fll8|) (with E\ = and 5 = 2e) were derived in |Q by solving equations 



for the superpotential which appear in the SUSY approach. In our terms, this approach 



deals with the function £(x) such, that £ has only one zero (just in the point Xo), £' changes 
sign nowhere and £ does not have poles on a real axis (otherwise they would give rise to 
additional zeros of W+). Therefore, in the situation considered in [|12|], ]TJ| corresponds to 
the ground state and i\) 2 describes the first excited state - in agreement with the conclusion 
of the previous section of the present article. Thus, in this particular case our approach 



reproduces the results of fll2| , |13 



V. ILLUSTRATIONS. DEFORMATIONS OF POTENTIAL LEAVING TWO 

LEVELS FIXED 

To illustrate the general results (0) - @, let us consider the following example: £ = x 4 + 
2x 2 x\ — x\. The derivative £' = at x = 0; therefore, as is explained in the preceding section, 
the corresponding example cannot belong to the set considered in ||12|| . After straightforward 
calculations, one obtains 

^ X 2 Xn 1 r . A\ Ao / \ 

U = E l + — -f + —AA a + 1 + 2 — , 20 

Ax\ Ax 4 x 2 + Xq (x 2 + Xq) 2 

where 

A = 2x 2 (2 + 4), At = (3xt + xi){5 - 4), A 2 = -(3x 4 + x 4 )x 2 (7 + 4)- (21) 

30-^ CC-^ X-^ 

The functions are equal to 

in = (x 2 + ^o)- Q exp(-^), a = (22) 



ip 2 = V>i (x 2 - x 2 ,) (a; 2 + a;^) , a; ± = ^/xq + xf ± a; 2 ,. (23) 



It is seen from (f22"|), ( |23|) that ^ has no nodes at the real axis, while ip 2 turns into zero 
at x = ±x_. Therefore, ipi corresponds to the ground state, while ip 2 describes the second 
excited state. 

As we see from ([5]), the Schwarzian derivative is an essential ingredient of the expression 
for the potential under discussion. It is known that the Schwarzian derivative is invariant 
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with respect to the linear-fractional transformations. Therefore, it is instructive to apply 
such a transformation to the potential as a whole and look at the resulting expression. Let 
us make the substitution 

£ = C2V + d2 . (24) 

CiT] + di 

We will use it below for generating in an explicit form rather rich families of the potentials, 
corresponding to two known levels. As \rf[ x remains invariant, only the part of (|8[) contains 
the terms proportional to AE and (AE) 2 , changes under this transformation. Then the 
potential and wave functions of the states under discussion take the form 

U = El -^ + ^ + ^ + ^-^Y -![,], (25) 
2 c\d 2 — c 2 di 4 rf l vf 1 2 

Y = AE (cir] + d 1 )( y c 2 r] + d 2 ) 

c\d 2 — c 2 d\ 

ipi,2 = e" p $i, 2 , $1,2 = c h2 7] + d h2 , (26) 
ji" — y 

P = 2 ^ , X = P - ln(ci?? + rfi)- 

If c 2 = = di, one can see that Y = AEr] and U(Ei, E 2 ,£) = U(Ei, E 2 , r/' 1 ) = 
U(E 2 , Ei,rj) in accordance with ([T4]). 
In the limit 

ci = = d 2 (27) 

we obtain the original potential U(Ei, E 2 ,£) = U (Ei, E 2l rf). 

Let us assume first we have some function r](x) characterized by the set of numbers 
(ni, n 2 , m* - )) introduced in Sec. III. The original potential has the form (|7|) with £ = 77. 



Then, let us take £(x) according to (33) with nonzero arbitrary coefficients q and dj. As the 



result of the transformation of (|24"D, each of the aforementioned numbers can change (for 
example, zeros scjj. of the combination C\T) + d\ generate poles of £, zeros Xj- 2 " 1 of c 2 r] + d 2 
correspond to zeros of £, each zero ^ of 7/ generates a zero of £'). Therefore, 
the levels £?2 which corresponded to the iVi-th and A^-th levels now can, in principle, 



correspond to another quantum numbers (Mi, M2). Making the transformation, inverse to 
(p4]), one may restore the form of the potential (|7]), but now £ 7^ r/, with £ having the form 
([24]), in which coefficients under discussion play the role of parameters. Thus, we obtain a 
family of deformations leaving two energy levels Ex, E2 fixed. These deformations can be 
described, on equal footing, by the deformation of the form of the function or of that 
of the potential. 

For definiteness, we will choose the second possibility. If ci, C2 7^ 0, one can always 
achieve c\ = c% = c by proper rescaling the function £(x) that does not affect, according to 



(jTjj), the function U(x). Then, defining 



c = 2(3, dx = -6- AE, d 2 = -5 + AE, 7 = (AE ]' - , (28) 

we obtain 

Y = (3n 2 - 8rj - 7, (29) 



V = El + M _ S + 2,3, - *[„]. + _ 4 W _ J, - 7 ), ,:«, 

2 2 4 77 2 77 



Below we will see how introducing nonzero parameters /3,7 affects the potential and wave 
functions (|26|). 



A. Example 1 

Let us choose 77 as a polynomial : 

77 = Aax(x — x ), a > 0, 77 = a(Vo + x — 2; ), 2; = x — x , Vq = const. (31) 
Demanding that p be regular at x = and at 2 = ±xo, we obtain from fl26|) the constraints 
7 = /V(K) + x A f + 8azjj - 6a(V Q + x 4 ) = /3aV 2 - Aax 2 - 5aV , (32) 

whence 
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v °-~^i~^h + w (33) 

7 = r\R ~j),R= ^ = ^ + 2« + 36x " 4 . 
The expression for the function £ reads 

2AE , . 

t = 1 + ^ + ^ + v (34) 



The potential can be obtained form (50) or directly from|7|. It has the form U = J2n=o c 2n% 2n , 
cio = ^|^-.It is convenient to rescale the variable in such a way that the coefficient in 
the potential U at the largest power be equal to 1. This can be achieved by x = Xy, 
(3aX Q = 8u, where to = 1 or — 1. After some manipulations we get the new potential 
U = X 2 U, corresponding to levels Ei t2 = X 2 E^ 2 : 

U = y 10 - Qfiy 8 + (13/x 2 + — )y & - (12/x 3 + 22u;)y 4 + (V + 31/zu; + -^)y 2 (35) 
, E 1 + E 2 15 2 

3fiuy 4 2 2 3 
P= ¥ —+y(»u + -), 

(AE) 2 = 16/i~ 2 (4/i 6 + 4,;/i 3 + 9). (36) 

where /i = x^/X 2 . (In fact, the formula ( |3~5| ) can be extended to negative /i as well.) The 
quantities /i and E 12 are not independent but connected by eq. (|36|) that appeared due to 
the condition of the regularity of the potential (|l^). It is worth noting that the parameter 
5 cancels and does not enter the expression (]35|) due to the conditions (|32|), (|33|) . 

It follows from (|26|), (f28|) and (|33"D that, up to the constant factor, the function $ 12 = 



z 2 -2{iz + qi j2 = {z-z 1 )(z-z 2 ), z = y 2 , z lj2 = fi±^Jfi 2 - q h2 , q h2 = \ + ± u. Let, 
for definiteness, E 2 > E\. Consider first the case uo = 1. Then after some algebra one easily 
finds that < q 2 < fi 2 and q\ > fi 2 . Therefore, the function $x does not have the nodes 
at the real axis and corresponds to the ground state. The function <3> 2 has four zeros and 
corresponds to the fourth state. In a similar way, we obtain that for uj = — 1 the quantities 
qi < 0, < q 2 < fi 2 , so the wave functions under discussion describe the second and fourth 
excited states. 
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B. Example 2 



One may exploit the ansatz (|TTD with £ = 77 for the potential (|30|) with ^7/0. 
Substituting it into fl26|), one obtains: 



' = "2 [L %^ " -T + V exp( V ^ - V exp( "V ^ )] " (37) 

Let us consider the example: 

W+ = ax + bx 3 , 7 = 0, 0^0, a = 5 = 5>0, b > 0. (38) 
After simple but cumbersome manipulations we get the potential 



U=^(- X j + 3r+ a ^-^)+u,r = fi^ t a* = °, (39) 
b 2 , ,ab P 2 , 4 (a 2 - 12b) 2 a 3 3 x 2 

U = —X 6 + (— + --)x 4 + - ; V-- + — -5 ^- + T- 



4 v 2 4 y 4 2 4(x 2 1 + x 2 ) 4(a; 2 + x§) 2 ' 

and wave functions 

, t = C+gT^ e-d - §,), * = , = 7 2L^, (40) 
/x 2 + xg)^P a ( x + ^x 2 + xl)^ ^Jx 2 + x 2 o 



[X 

It is seen from ( [4"0"| ) that rj' > and the potential (|3Up is regular for any choice of parame- 
ters, so the conditions fllED , (|i~6[) are irrelevant for this case. The functions ipi and ip2 are 
normalizable, provided (3 2 < ab. It can be readily seen from ( |40f) that the function ip 2 has 
one node at x — 0, whereas ^1 turns into zero nowhere. Thus, ip\ and ip2 correspond to the 
ground and the first excited states, respectively. In the limit (3 = we reproduce the result 



for the example 3 of [12 . 



C. Example 3 



Let now W+ = A(shx — shx ), 7 = 0, (3 ^ 0, 5 = 5. Then, repeating calculations for 
this case, we get 

12 



2 2 c hx ch^^ 2chx ch 2 x 

^o(^) = "7X5 — ( shx ~ shx ) 2 - — — (2chx - chx ) + -, 

a = (2c/ix ) _1 [5c/?.x — j3sh(x — x ) + x{(5 — Sshx )]. 

Here ?7o is the potential corresponding to the anisotropic paramagnet of the spin 1/2 in an 
oblique magnetic field 0. The function rf > 0, so p(x) is regular in any point for any choice 
of parameters. The wave function is normalizable provided — 5e~ x ° < j3 < 5e x °. One can 
easily show that it follows from this condition that ipi does not have nodes and corresponds 
to the ground state, while ip2 has one node at x = xq and corresponds to the first excited 



state. In the limit f} — the example 1 of [12] is reproduced 



VI. CONCLUDING REMARKS 

Thus, in a very simple and direct approach we found a rather general solution that gives 
us the structure of potentials with two known eigenstates E\, E2 in terms of one function £(x) 
and one parameter coinciding with the energy difference AE — E 2 — E\. Moreover, we get 
not only the potential itself but, also (in terms of the SUSY language), the superpotential. 
Depending on properties of the function £ (x) and the type of the regularity condition of the 
potential in the vicinity of zeros of £'(x) ( dl5|) or (|16|)), one can obtain not only the ground 
or first excited state but, in principle, any pair of levels. The natural question arises whether 
the approach of the present paper is extendable to the case of three (or more) levels. This 
problem needs separate treatment. We hope that movement in this direction will promote 
further understanding links between QES-type systems, SUSY and the inverse scattering 
approaches. 
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acknowledges gratefully finansial support from the German Academic Exchange Service 
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